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ABSTRACT 

Material ejected during (or immediately following) the merger of two neutron stars may assemble 
into heavy elements by the r-process. The subsequent radioactive decay of the nuclei can power 
electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying 
such events is an important goal of future transient surveys, offering new perspectives on the origin 
of r-process nuclei and the astrophysical sources of gravitational waves. Predictions of the transient 
light curves and spectra, however, have suffered from the uncertain optical properties of heavy ions. 
Here we consider the opacity of expanding r-process material and argue that it is dominated by line 
transitions from those ions with the most complex valence electron structure, namely the lanthanides. 
For a few representative ions, we run atomic structure models to calculate radiative data for tens 
of millions of lines. We find that the resulting r-process opacities are orders of magnitude larger 
than that of ordinary (e.g., iron-rich) supernova ejecta. Radiative transport calculations using these 
new opacities indicate that the transient emission should be dimmer and redder than previously 
thought. The spectra appear pseudo-blackbody, with broad absorption features, and peak in the 
infrared (~ 1 (an). We discuss uncertainties in the opacities and attempt to quantify their impact on 
the spectral predictions. The results have important implications for observational strategies to find 
and study the radioactively powered electromagnetic counterparts to compact object mergers. 



1. INTRODUCTION 



Hydrodynamical simulations suggest that a small frac- 
tion of mass is ejected when two neutro n stars (or a black 
hole and neutron star) collide or merge flJanka et al.|1999| 



Rosswog ct al. 1999; Lee 2001; Rosswog 2005; Occhslin 
et aL | 20071|Chawla et al.|2010||Shibata fc T aniguchi|2011| 
Hotokezaka et al.||2013p . If this ejecta is sufficiently neu- 
tron rich, it will assemble within seconds into hea vy el 



timer fe Schramm|19"74||Eichler et aL|1989||Freiburghaus 



ements via rapid neutron captures (the r-process) (|Lat 

Jiliii 



et aL 1999p . The subsequent beta decay of the nu- 
clei will heat the ejecta for days, powering a thermal, 
supernova- like transient ( Li fc Paczyhski|1998 ). Because 
the ejected mass is small in comparison to ordinary su- 
pernovae (SNe), the light curves of these "r-process su- 
pernovae" are expected to be relatively dim and short 
lived. Previous radiative models predict peak bolomctric 
luminosities around 10 40 — 10 42 ergs s _1 , peaking at opti- 



cal wavelengths and lasting around a day (|Li fc Paczyhski 



1998 



2011 



Kulkarni 2005; Mctzger et al 



Goriely et al.||2011| |Piran et a 



20101 IRoberts et al. 



2012) 



XTthougfi we have not yet discovered an r-process SN 
from a neutron star merger (NSM) , there are compelling 
reasons to look for them. Because these outflows are 
non-relativistic, they emit radiation relatively isotropi- 
cally, and are therefore promising electromagnetic coun- 
terparts to gravitational wave sources; if discovered coin- 
cidently, they could enhance the scientific value o f an ad- 
vanced LIGO/VIRGO gravitat i onal wave signal ([S chutz 
19861 IKochanek fc PiranTT993l ISylvestre|[2003l iPhinney 
2009 Mandel fc O'Shaughnessy|20101 |Metzger fc Berger] 
2012||Kelley et al.|2012||Nissanke et al.|2012[ ). Discovery 
ot r-process SJNe would also dramatically illuminate our 
incomplete understanding of heavy element production 



in the Universe. The NSM ejecta is thought to be a re- 
markably pure sample of r-process material, which would 
allow us to cleanly study heavy elements near their pro- 
duction site, and soon after they had been created. In 
principle, analysis of the light curves and spectra of these 
radioactive transients could be used to quantify the mass 
and chemical composition of the ejecta, which would clar- 
ify the unknown site(s) of r-process nucl eosynthesis (e.g., 
Arnould et al.||2007| [Sneden et al.||2008[ ). 



Perhaps the largest remaining uncertainty in our un- 
derstanding of r-process SNe has concerned the opacity 
of the ejected debris, which (along with the ejecta mass 
and kinetic energy) is a key parameter determining the 
brightness, duration, and color of the transient. The 
ejecta of NSMs consists of heavy elements in rapid differ- 
ential expansion, and at relatively low densities and tem- 
peratures (p ~ 10~ 13 g cm~ 3 and T ~ 5000 K at 1 day 
after ejection). Because almost nothing is known about 
the optical properties of such material, previous radiative 
transfer models have simply adopted opacities character- 
istic of ordinary SNe. In Type la supernovae (SNe la), 
for example, the opacity is due primarily to numerous 
iron group lines, which are blended by Doppler broaden- 
ing into a pseudo-continuum. We can expect that lines 
will also dominate the opacity of NSM ejecta, but un- 
fortunately very little atomic data exists for ions heavier 
than the iron group, either from theory or experiment. 

Given the ignorance, we might first consider some gen- 
eral expectations from atomic physics. The number of 
strong lines will be larger for ions with greater complex- 
ity - i.e., with a denser packing of low lying energy states. 
Naively, one might expect higher Z elements to be more 
complex than the iron group. Of course, what matters 
is not the total number of electrons, but the number of 



2 



Kasen, Badnell, and Barnes 




40 60 

atomic number Z 

Fig. 1. — Complexity of the elements (top panel) and their mass 
fractions in the r-process ejecta of neutron star mergers (bottom 
panel). The top panel plots the number of states in the ground 
configuration for singly ionized ions, as estimated using the simple 
permutation counting of eq. ^ The pattern of peaks reflects the 
filling of valance shells, with the color shading giving the orbital 
angular momentum I (yellow = s, blue = p, green = d, red = /). 
The bottom panel plot s the mass fractions o f the tidal tail ejecta 
from the simulations of |Roberts et ah] ( |201 1 [ > . 

distinct ways of distributing valence electrons within the 
open shells. A subshell with orbital angular momentum 
I has g = 2(2Z + 1) magnetic sublevels; one can estimate 
the number of states in a given electron configuration 
by simply counting the permutations of the valence elec- 
trons 



c = n, 



rii\(gi - Ui)V 



(1) 



where n, is the number of electrons in the n^-orbital la- 
beled i, and the product runs over all open shells in 
a given configuration. The different terms and levels 
(i.e., distinct combinations of L, S, J) derived from these 
various permutations are split by electrostatic and fine- 
structure (e.g. spin-orbit) interactions. Equation Jl] can 
thus be used to estimate the relative number of distinct 
energy levels of an ion, while the number of lines (i.e., 
radiative transitions between levels) will scale roughly as 
C 2 . Figure [l] plots the complexity measure C for the 
ground configurations of singly ionized ions, where the 
pattern of I shell filling is clearly seen. 

Equation [T] provides immediate insight into the opac- 
ity of r-process ejecta. Ions with valence shells of higher 
/ are more complex, as are those whose open shells are 
closer to half filled. This is why the iron group, with 
a nearly half-filled d (I — 2) shell, usually dominates 
the line opacity in typical astrophysical mixtures. Heavy 
r-proccss ejecta, however, includes uncommon species 
of even greater complexity. Of particular importance 
are the lanthanides (58 < Z < 70) and the actinides 
(90 < Z < 100) which, due to the presence of an open / 
(I = 3) shell, have complexity measures roughly an order 
of magnitude greater than the iron group. While the ac- 
tinide series is generally of very low abundance, the lan- 
thanides may represent several percent of r-process ma- 
terial by mass. We will find that these species dominate 
the total opacity of NSM ejecta, resulting in opacities 
orders of magnitudes greater than previously assumed. 



To calculate the ejecta opacity in detail, we need a 
comprehensive list of atomic lines. As almost no data is 
available for heavy ions, we turn here to ah initio a tomic 
structure modeling using the Autostructure code (Bad- 
nell|2011 ) . These models determine the approximate ion 
energy level structure and the wavelengths and oscilla- 
tor strengths of all permitted radiative dipole transitions 
(fj3|. Without fine tuning the structure model, the com- 
puted energies and line wavelengths are not exact. For- 
tunately, the effective opacity in an expanding medium 
is a wavelength average over many lines. Because our 
models reasonably capture the statistical distribution of 
levels and lines, they can be used to derive fairly robust 
estimates of the pseudo-continuum opacity (q3q5|). 

Modeling the radiative properties of all hign Zions is 
a long term endeavor; we present here initial structure 
calculations for a few representative ions selected from 
the iron group (Fe, Co, Ni), the lanthanides (Ce, Nd), 
and a few other heavy d-shell and p-shell ions (Os, Sn). 
The Autostructure line data is then used to calculate 
the opacity of expanding ejecta under the assumption of 
local thermodynamic equilibrium (LTE). We show that 
ions of similar complexity have similar properties, which 
allows us to estimate the total opacity of an r-process 
mixture based on the representative species (^j6|. 

The derived opacities can be input into a multi- 
wavelength, time-dependent radiative transfer code to 
predict the observable properties of r-process supernovae 
(Q. We discuss here the general spectroscopic proper- 
ties of these transients, while a companion study explores 
the broadband lig ht curves and their de pendence on the 
ejecta properties ( jBarnes k. Kasen|2013 |. In general, the 
high r-process opacities result in light curves that are sig- 
nificantly broader, dimmer, and redder than previously 
believed. These results have important implications for 
observational strategies to find and interpret the radioac- 
tively powered electromagnetic counterparts to NSMs. 

2. OPACITY OF RAPIDLY EXPANDING EJECTA 

We set the stage by reviewing the physical properties 
of the material expected to be ejected in NSMs. We 
then describe the nature of the opacity in such gas, in 
particular that arising from line interactions in a rapidly 
expanding medium. 

2.1. Physical Conditions of the Ejecta 

There are at least two distinct mechanisms by which 
material may be ejected in NSMs: 1) During the merger, 
surface layers may be tidally stripped and dynamically 
flung out in "tidal tails" . 2) Following the merger, mate- 
rial which has accumulated in a centrifugally supported 
disk may be blown off in a neutrino or nuclear driven 
wind (Lcvinson 2006} [Surman et aL]|2006 Metzger et al. 
|2008||20U9[ ). The amount ot mass ejected in the tidal tails 
appears to depends upon many factors: the NS mass ra- 
tio, the equation of state of nuclear matter, and the treat- 
ment of gravity, but simulations give values in the range 
M e j = 10~ 4 - 10 _1 M & . A similar amount of mass may 
potentially be ejected in the disk wind. In both cases, 
the characteristic velocities are w C j ~ 0.1 — 0.3c, of order 
the escape velocity from a NS. 

The composition of the material ejected by the two 
mechanisms is likely different. The tidal tail ejecta is 
initially cold and very neutron rich (electron fraction 
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Y e ~ 0.1), and should rapidly assemble into heavy el- 
ements (Z > 50) through the r-process. The conditions 
in the disk wind are quite different; weak interactions 
will drive the material to be less neutron rich (Y e ps 0.5) 
and the entropy will be higher. This environment is more 
similar to the neutrino driven wind from proto-neutron 
stars in core collapse supernovae. It is unclear whether 
a robust r-process occurs in such a wind, or whether the 
distribution only extends to atomic numbers Z ~ 50. If 
neutrinos drive Y e close to 0.5, the composition ma y be 
dominated by radioactive 56 Ni ( Surman et al.|[2008 ). 

Soon after the mass ejection (~ 100's ot seconds) hy- 
drodynamical and nucleosynthetic processes abate, and 
the ejecta reaches a phase of free-expansion. In the ab- 
sence of any forces, the velocity structure becomes ho- 
mologous - i.e., the velocity of any mass element is pro- 
portional to radius, v = rt e j, where t e j is the time since 
homology sent in. The density of the ejecta at 1 day can 
be estimated 
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2.8 x 10 
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g cm \ 



(2) 



where M_ 



A/ t ,j/lll 2 M, 



O, wo.i = v ei /0.1c, and t d 



i e j/ day. In just a day, the density of the ejecta has 
dropped by ~ 20 orders magnitude from its original value 
in the neutron star. 

The ejecta material initially cools very effectively by 
expansion, but will be reheated by the decay of r-process 
nuclei. Radioactive energy is released in the form of 
gamma-rays, beta particles, and fission fragments, which 
will be thermalize d, to various degre es, by scattering 



within the ejecta ( Metzger et al.[2010[ ). The heated ma- 
terial will radiate, and thermal photons will escape the 
medium on the effective diff usion timesca le for a homol- 
ogously expanding medium ( Arnett||l980 ) 



M cj K 
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1.7 M 
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where the opacity k of the ejecta has been normalized 
to a value Kq.i = ft/0.1 cm 2 g _1 , a value appropriate for 
iron group elements (but not, we will find, for r-process 
elements). This timescale for diffusion sets the duration 
of the radioactively powered light curve. 

The luminosity near the peak of the light curve will 
be of order the instantaneous rate of energy deposition 
L ps M e je, where e is the r adioactive ene rgy released per 
unit time per unit gram (Arnett 1982). Stefan's law, 
L = 47rr 2 (7 s bT 4 , provides an estimate of the surface tem- 
perature 
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V&MH^voxtd)- 1 ' 2 K. (4) 



For times t cj > 1 day, the characteristic temperatures 
and densities of NSM ejecta are thus roughly in the range, 
T ~ 10 3 — 10 4 K and p ~ 10~ 16 - 10~ 12 g cm" 3 . Un- 
der these conditions, and assuming local thermodynamic 
equilibrium (LTE), the gas will be in a low ionization 
state, being mostly singly or doubly ionized near the pho- 
tosphere. 

2.2. Local Thermodynamic Equilibrium 



We adopt LTE to compute level populations in this 
paper, a necessary approximation given the complexity 
of the ions involved. The low density of NSM ejecta 
a t t e j > 1 day is not sufficient for collisional processes 
alone to establish LTE. However, in the optically thick 
regions below the photosphere, the radiation field will 
tend toward a blackbody distribution and radiative tran- 
sitions will drive the level populations to their LTE val- 
ues. Because the effective diffusion time and spectral 
energy distribution (SED) are mostly set by processes 
near and below the photosphere, LTE calculations likely 
provide a reasonable first approximation. At late times 
(t > 20 days), when the entire remnant becomes trans- 
parent, LTE will break down at all radii and result in 
poor SED predictions. 

Considering the heavily radioactive environment of 
NSM ejecta, one may worry that, even at early times, de- 
partures from LTE may be driven by non-thermal ioniza- 
tion/excitation processes (namely, impacts by fast elec- 
trons that have been Compton scattered by radioactive 
gamma-rays). As a rough estimate of the potential ef- 
fects, we compare the rates for a bound-bound transi- 
tion of energy AE. The non-thermal excitation rate is 
R nt ps fe/AE, where e is the radioactive power released 
per particle, and / is the fraction of that power that goes 
into excitation (as opposed to ionization or thermaliza- 
tion). The radiative excitation rate, assuming a black- 
body field, is Rbb = B\2By{T) where B12 is the Einstein 
absorption coefficient and B V (T) the Planck function. 
Using the Einstein relations, the ratio of rates is 



-Rbb 



AEA 21 



(e AE / kT -l), 



(5) 



where the Einstein spontaneous emission coefficient 
A21 ~ 10 8 — 10 9 for permi tted optical tran s itions . For 



r-process ejecta at 1 day, |Metzger et al. (2010) find 
e ps 1 eV s" 1 Considering values AE ~ 1 eV, T ps 5000 K 
and / ~ 1/3, the ratio -R n t/-Rbb ~ 10 -8 . Non-thermal 
transitions are therefore negligible except for transitions 
well above the thermal energy (AE/kT > 20). A sim- 
ilar argument can be made for non-thermal ionization. 
We conclude that the radioactive energy deposition will 
likely not seriously undermine our LTE assumption, ex- 
cept at late times when the ejecta becomes rather cold 
and transparent. 

2.3. Line Expansion Opacity 

The opacity of bound-bound transitions is significantly 
enhanced by the high expansion velociti es found in super- 
nova and NSM ejecta ( Karp et al.|1977 ). We discuss here 
the "expansion opacity" formalism which wavelength av- 
erages the contribution of multiple lines and treats the 
line radiative transfer in the Sobolev approximation. 

As a photon propagates through the differentially ex- 
panding medium, its wavelength is constantly Dopplcr 
shifted with respect to the comoving frame. For a ho- 
mologous (Hubble-like) expansion, this Doppler shift is 
always to the red, and proportional to the distance trav- 
eled. The photon will interact with a line when its co- 
moving frame wavelength is redshifted into resonance 
with the line rest wavelength. The spatial extant of this 
interaction, or resonance, region is set by the intrinsic 
width of the line profile. If, for example, the line width 
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is due to the thermal velocity, Avt, of ions, the physical 
size of the resonance region is As ~ Avtt e - 3 . 

Because the thermal velocities in NSM ejecta are very 
small (v t ~ 1 km s _1 ) compared to the ejecta velocities 
(i> e j PS 10 5 km s _1 ) the resonance region of a line is tiny 
compared to the ejecta scale height, and the matter prop- 
erties can be taken to be constant over the region. This i s 
the essence of the Sobolev approximation (Sobolev 1960). 
In this limit, the line extinction coefficient can be analyt- 
ically integrated to give the Sobolev optical depth across 
the resonance region 



T s — /osc^l^cj Ao, 

m e c 



(6) 



where / osc is the oscillator strength and Ao the rest wave- 
length of the line. The Sobolev optical depth is a local 
quantity which depends on n\, the number density in the 
lower level of the transition at the location of resonance. 
The probability that a photon interacts (i.e., is scattered 
or absorbed) a least once in traversing the resonance re- 
gion is simply 1 — e~ T " . 

A photon traveling through the expanding medium 
comes into resonance with lines one-by-one, sweeping 
from blue to red. The effective mean free path depends 
not on the strength of any one line, but rather on the 
wavelength spacing of strong (r s » 1) lines, which can 
be quantified as follows. Say that within some wave- 
length region (A,, Aj + AA<) we have N strong lines. The 
spacing between the lines is, on average, AAj/iV and the 
velocity gradient which Doppler shifts a photon from one 
line to the next is Av s /c = AXi/XiN. If homologous ex- 
pansion holds, the distance a photon travels between line 
interactions is then l m f p = Av s t c j. This is an estimate 
of the mean free path, while the inverse quantity, i~f p , 
defines the matter extinction coefficient (units cm -1 ) 
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where the bin size AA^ can be chosen arbitrarily to av- 
erage over a reasonable number of lines. 

A formal deriva tion along these lines was introduced by 
Karp et al. ( 1977 ) to estimate the extinction coefficient in 



an expanding medium. We use here the slightly m odified 
expression developed by |Eastman fc Pinto| ( |1993"| ) 



(S) 



where the sum runs over all lines in the wavelength bin 
AAi . Eq. M takes into account the cumulative effect of 
many weatc lines; in the case that all lines are optically 
thick (r s 3> 1), it reduces to the simple estimate eq. [7 

An interesting property of the expansion opacity is the 
dependence on ion density, which appears only in the 
Sobolev optical depth, r s . In the limit that all lines are 
weak, the density dependence is linear. However, in the 
opposite limit where all lines are very optically thick, the 
extinction coefficient is independent of density. In our 
actual calculations, the density dependence (which is set 
by the statistical distribution of line optical depths) is 
weak, roughly logarithmic. As a consequence, certain 
ions can make a significant contribution to the opacity 
even when their abundance is very low. 



The expansion opacity (units cm 2 g x ) of bound- 
bound transitions is related to a ox by 



n cx (X) 
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Given the weak dependence of a cx on density, the line 
expansion opacity actually increases as the density drops. 

2.4. Applicability of the Sobolev Approximation 

At least three conditions must be met for the Sobolev 
approximation, which underlies the expansion opacity 
expression eq. [9j to be valid. The first, already men- 
tioned, is that the thermal velocity of the ions (pre- 
sumed to set the intrinsic line widths) must be signif- 
icantly smaller than the velocity scales over which the 
ejecta properties vary. For NSM ejecta, the ratio is 
Vt/v c j f=a 10~ J <C 1, which assures the applicability of 
the narrow line limit. 

A second condition, relevant for time- varying environ- 
ments, is that the time photons spend within a resonance 
region be short compared to the timescale over which 
the ejecta properties vary (e.g., the expansion timescale). 
For strong lines, a photon may scatter multiple times 
within the resonance region before finally being redis- 
tributed to the edge of the line profile and escaping (i.e., 
redshifting past) the transition. In the Sobolev formal- 
ism, the probability that a photon escapes the line after 
a scatter is 



(3 = 



1-e 



r c 1 for t > 1, 



(10) 



and the average number of scatters in a thick line is N sc ~ 
1//3. Assuming the distance traveled between scatters is, 
on average, wt^cj, the time spent in the resonance region 
is t sc ~ N sc vtt c j/c. The condition t sc < i e j sets a limit 
on the optical depth of the line r s < c/v% = 3 x 10 5 t! t x, 
where v t .i — v t /l km s _1 . In practice, optical depths 
> 3 x 10 5 are regularly encountered, in particular for 
resonance lines. In most cases, however, fluorescence 
provides an avenue for escape. The probability of de- 
excitation to a lower level is suppressed by a factor of /?, 
and so it is likely that the ion will eventually de-excite 
via a cascade through multiple low r s transitions. This 
generally evades the problem of extended line trapping, 
except perhaps for those few transitions in which fluores- 
cence is not possible. 

A third condition is that the intrinsic profiles of strong 
lines must not, in general, overlap, as this would intro- 
duce a coupling of the radiative transport between lines. 
In particular, a photon that escapes from the red edge of 
one line will have an enhanced probability of escaping a 
second overlapping line. This invalidates the sum in eq.|9] 
which assumes an independent interaction probability for 
each line. Overlap of weak lines (which are extremely 
numerous in our calculations) is common, however this 
likely does not introduce any serious error, as the t s de- 
pendence is linear when r s < 1. Moreover, the opacity is 
usually dominated by the strong lines. Occasional over- 
lap of strong lines is inevitable, and may moderately re- 
duce the expansion opacity at certain wavelengths. The 
entire Sobolev formalism, however, becomes inapplicable 
when the wavelength spacing of strong lines, AX/N ', be- 
comes comparable to the intrinsic (e.g., thermal) width 



r-process opacities 



■5 



AAt = Ao(i>t/c) of the lines. From eq. |9]we can define 
opacity when such "saturation" occurs 
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(11) 

When k ox > K sat , strong line overlap is the norm and 
the Sobolev expansion opacity formalism can no longer 
be trusted to return reasonable values. Under some con- 
ditions, and at certain wavelengths, we will find that our 
calculated r-process opacities approach or exceed satura- 
tion, such that this issue may be a serious concern. 

2.5. Other sources of Opacity 

Other potential sources of opacity include free-free 
(i.e., bremsstrahlung) , bound-free (i.e., photoionization) , 
and electron scattering. In NSM ejecta, none of these 
turn out to be important compared to bound-bound. For 
example, the wavelength independent electron scattering 
opacity is given by 

« cs = ^0.4(f)cm 2 g- (12) 

where x is the mean ionization fraction and A the mean 
atomic weight of the ions. For NSM ejecta comprised of 
lowly ionized (x ~ 1) heavy elements (A ~ 130), k cs is 
a factor x/A ~ 10~ 2 smaller than the typical value for 
ionized hydrogen, and much less than the r-process line 
opacity at all wavelengths of interest. 

The free-free opacity for a gas i n ionization state x 
1 is g iven approximately by (e.g 
19861 



Rybicki & Lightman 



Kff - 0.005^p_ 13 T 4 1/2 (A/l M m) 



3 cm 2 g 1 , 



(13) 



where we have set the correction for stimulated emis- 
sion and the Gaunt factor to unity. For the low den- 
sities and x/A values found in NSM ejecta, the free- 
free opacity at the relevant wavelengths is negligible 
(ka w 1CT 6 cm 2 g" 1 ). 

Finally, the opacity due to a bound-free transition for 
photons at the threshold energy (where the cross-section 
is largest) for some excited level of an ion is given by 



ttbf 



o-o e - AE / kT 
Am p Z(T) 



(14) 



where AE is the excitation energy of the level and Z(T) 
the partition function (LTE is assumed) . In order for the 
bound- free transition to apply to optical/infrared pho- 
tons, the level must be a highly excited state, AE > 
X — 2 eV, where x is the ionization potential. As a 
representative estimate of the opacity at T = 5000 K, 
we consider a singly ionized heavy element with A = 
120, AE ~ 8 eV and Z(T) = 20 and adopt the hydro- 
genic value <7n ~ 6 x 10~ 18 cm 2 . The resulting opacity at 
threshold is also very small, «hf ~ 1*3 x 10~ 5 cm 2 g _1 , 
due mainly to the Boltzmann factor. At ultraviolet wave- 
lengths (A < 1000 A) the bound-free opacity may actu- 
ally dominate, since for low lying levels (AE ps 0) one 
finds Kbf 
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Fig. 2. — Atomic structure model calculations of the excitation 
energy of the lowest level of Ndll electron configurations. Cir- 
cles denote the results from Autostructure obtained under various 
optimization approaches (described in text). Stars denote the ex- 
perimental energies from NIST. 



3. ATOMIC STRUCTURE CALCULATIONS 

To estimate the radiative data fo r high-Z elemen ts, we 
used the Autostructure program (Badnell 2011 1. This 
code was used previously to calculate data up t o Ni for 
the update d opacities of the Opac ity Project (|Badnell| 
et al. 2005| ). Recent developments ( Badnell et al.|[2012| ) 
have enabled it to be used to make extensive calcula- 
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tions of radiative (and autoionization) rates for a half- 
open f-shell. Autostructure calculates the approximate 
level energy structure of ions, and all relevant radiative 
transition rates, given a user specified set of electron 
configurations. The many-electron quantum mechani- 
cal problem is treated using a multi-configuration wave- 
function expansion with a Breit-Pauli Hamiltonian. We 
used the kappa-av eraged relativistic wavefu nction option 
as introduced by |Cowan fc Griffin (1976). The radial 
orbitals were determined using a Thomas-Fermi-Dirac- 
Amaldi potential. As the standard LS'-coupling scheme 
breaks down for high Z elements, we adopted a level- 
resolved intermediate co upling scheme. 

We used the NIST ( |Kramida et al.| |2012[ ) atomic 
database to identify the electron configurations corre- 
sponding to the ground and low-lying states of each ion. 
For several of the high-Z ions, the NIST data appeared to 
be incomplete, and we included additional configurations 
suspected to be relevant. Table [T] lists the configurations 
used for each ion. For the lanthanides, the highest nl 
orbital we considered was the 6p one. We experimented 
with including configurations generated by electron pro- 
motion to higher orbitals (e.g., n — 7,8), however these 
typically produced highly excited levels not significantly 
thermally populated under the relevant physical condi- 
tions. As we did not did not notice any large effects on 
the opacities, we omit these configurations from our fi- 
nal calculations, although more exhaustive explorations 
of configuration space are certainly warranted. 

Autostructure includes a dimensionless radial scaling 
parameter for each n^-orbital, which must be optimized 
to establish a realistic level structure for low-charge ions. 
The optimization consists of varying the scaling param- 
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TABLE 1 

AUTOSTRUCTURE ATOMIC STRUCTURE MODELS 
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eters so as to minimize a user specified weighted sum of 
eigenenergies. The closed-shell core cannot be excluded 
from the structural optimization of complex heavy near- 
neutral ions because of the str ength of core p olarization 
effects on the valence orbitals ( Palmeri et al.pOOO ). We 
therefore used a single common variational scaling pa- 
rameter for all closed-shell orbitals, but varied the pa- 
rameters of the valence orbitals independently. This ab 
initio optimization procedure does not require any ob- 
served energies. Thus, it is ideally suited to situations 
such as the present one where the observed data is at 
best sparse. 

We explored several strategies for optimization. The 
first - which we label optl - was to simultaneously 
vary the scaling parameters for all included (core-plus- 
valence) orbitals and to minimize the equally weighted 
sum of all energy levels included by the configuration 
expansion. This has the advantage of not biasing the 
structure towards any given configuration(s) , which is 
valuable given that we seek radiative data for many ex- 
cited levels. The disadvantage is that it gives no due 
preference to the ground state, and thus does not always 
predict a ground state configuration in agreement with 
what is deduced experimentally. Figure [2] shows the en- 
ergy level structure for Ndll, where it is seen that the 
optl optimization predicts the wrong ground configura- 
tion. While for many applications this would be a fatal 
flaw, in the present context we expect several low ly- 
ing configurations to be significantly populated, in which 
case the mean line opacity may be less sensitive to the 
exact configuration ordering. 

We considered a second strategy (opt2) whereby the 
above optimization was first applied to only those or- 
bitals included in the ground configuration. These scal- 
ing parameters were then fixed, and a second optimiza- 
tion was carried out varying the parameters of all remain- 
ing orbitals. This method usually produced the correct 



ground state configuration. The energies of the excited 
levels were also close to but a bit higher than the avail- 
able NIST values, and overall not as good as those found 
using the optl approach (Figurej2j). 

The model structure can be further refined by itera- 
tively adjusting the scaling parameters by hand. We at- 
tempted this for Ndll, guided by the trends found in the 
optl and opt2 calculations. An improved solution was 
found (optS) which reproduced the ground and first two 
excited level energies almost exactly. Further iterations 
could presumably improve the result, but this sort of 
manual alignment is time consuming, and more art than 
science. We only attempted this optS approach only for 
Ndll, which is the most important ion for our r-process 
light curve calculations. 

4. IRON GROUP OPACITIES 

4.1. Comparison to Kurucz Line Data 

The atomic properties of Z < 30 ions are reasonably 
well known based on experiment and previous struc- 
ture modeling. In particular, R. Kurucz has gener ated 
extensive lin e lists, including CD23, (~ 500, 000, |Ku 



Kurucz (1993 



rucz fc Bell (U995aJ), and CD1 (~ 42 million fines 



ese lists (which are dominated by 
iron group lines) have been derived from atomic struc- 
ture calculations using the Cowan code (Cowan 19811 
which have been iteratively tuned to reprod uce the ex- 
tensive obs erved experimental level energies ( |Kurucz fc 
Bell 1995b I. Supernova modelers have used the Kurucz 
data to successfully model the optical ligh t curves and 



spectra of observed (ir on-rich) SNe la (e.g., Kasen et al 



12009, 



Sim et al. 2010[ ) which suggests that, lor the iron 



group, the Kurucz line data can be taken to be reason 
ably accurate and complete. 

To validate our ab initio Autostructure line data 
against the observationally constrained data of Kurucz, 
we ran structure models for the first 4 ionization stages of 
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Fig. 3. — Calculated line expansion opacities for a mixture of iron group elements (10% Ni, 80% Co, 10% Fe) representative of decayed 
56 Ni. The right panel plots the Planck mean opacity versus temperature, while the left panel plots the wavelength dependence of the 
opacity at a specific temperature, T = 5000 K. The opacities calculated using our Autostructure derived line data generally agree with 



those using the Kurucz linelists. The calculations assume p = 10 g cm 



10 days, and a wavelength binning AA = 0.01A. 



Fe, Co, and Ni, using the electron configurations listed 
in table [I] Unlike Kurucz, we made no attempt (be- 
yond our ah initio optl optimization scheme) to tune the 
model, and our calculated level energies can differ from 
the experimental values by factors of 2 or more. Nev- 
ertheless, we find that our derived iron group expansion 
opacities are in good agreement with those of Kurucz 
(Figure |3|. Our Planck mean opacities differ from those 
of Kurucz by only ~ 30% over the temperature range 
1000 — 20, 000 K, and the wavelength dependence of the 
opacity is quantitatively similar, with the opacity rising 
sharply to the blue. The good agreement indicates that 
our Autostructure calculations capture the statistical 
properties of the lines, even if the individual energy lev- 
els and line wavelengths may be inaccurate. 

4.2. Application to Supernova Modeling 

To demonstrate how our Autostructure derived iron 
group opacities perform in a real transport calculation, 
we calculated synthetic light curves and spectra of a sim- 
ple SN la model. As a numerical "thought" experiment, 
we pretended that our understanding of SNe la was as 
rudimentary as it is for NSM outflows, and that our only 
expectation was that a carbon/oxygen star was burned 
to nuclear statistical equilibrium (NSE). We thus con- 
structed a spherically symmetric ejecta model consisting 
of uniform density, initially pure 56 Ni with a total mass 
of 1 Mq and kinetic energy of 10 51 ergs, roughly the 
nuclear energy released in burning the C/O to NSE. Of 
course, real SN la are not homogenous, and in addition 
to 56 Ni are observed to contain a significant amount of 
intermediate mass elements (IMEs, Si, S, Ca). 

We calculated synthetic observables of this ejecta 
model using the time depend ent, multi-wavelen gth radi- 
ation transport code Sedona ( Kasen et al.|2006 l, and as- 
suming that the level populations were in LIE. Figure [4] 
shows that, despite the simplistic nature of the ejecta 
model, the broadband light curves qualitatively resemble 
those of observed SNe la, peaking at a B-band magnitude 
around -19 about 20 days after explosion. Overall, the 
light curves calculated using our Autostructure derived 



linelist are rather similar to those calculated using the 
Kurucz linelist, although differences up to 1 magnitude 
are seen at some epochs. 

Figure ]4] shows that the model spectrum (at 50 days 
after explosion) derived from the Autostructure linelist 
is also similar to that using Kurucz. Both calcula- 
tions resemble the SED of an observed SN la. The 
Autostructure model does not reproduce the positions 
of most spectral features, which is to be expected given 
that the line wavelengths are only approximate. Even 
the Kurucz calculation fails to reproduce every observed 
spectral feature, as the underlying ejecta model did not 
include the IMEs present in real SNe la. 

These results indicate that line data derived from our 
Autostructure models can be used to predict supernova 
SEDs (but not line features) with some reliability. The 
general agreement of our synthetic observables with those 
of real SNe la suggests that - even with very crude knowl- 
edge of the underlying ejecta structure - we may still be 
able to predict the light curve and colors of radioactive 
transients to a reasonable level of accuracy. 

5. HIGH Z OPACITIES 

We have calculated structure models for several ele- 
ments beyond the iron group, including tin (Sn, Z = 50, 
p-shell), cerium (Ce, Z = 58, f-shell), neodymium (Nd, 
Z = 60, f-block), and osmium (Os, Z — 76, d-shell). 
These species were chosen to sample different blocks on 
the periodic table corresponding to valence shells of dif- 
ferent orbital angular momentum. The total number of 
atomic levels/lines determined by the structure models 
are listed in Table 1 and illustrated in Figure [5j and are 
generally consistent with the simple complexity estimates 

of© 

As expected from simple physical arguments, we find 
that more complex atoms, in particular the lanthanides, 
have higher line expansion opacities. Figure [6] shows 
that the Planck mean opacity of neodymium is a fac- 
tor ~ 10 — 100 greater than that of iron, depending on 
the temperature. This is roughly consistent with the 
estimate one gets by squaring the complexity measure 
(equation [T]) to gauge the relative number of strong lines, 
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Fig. 4. — Synthetic light curves and spectra of a model consisting of pure 56 Ni ejecta, computed with different line data. The left panel 
shows that the SED (at 50 days after merger) calculated using the Autostructure line data resembles that obtained using the Kurucz data. 
Both resemble the spectrum of the Type la SN 2003du observed 34 days after peak (~ 52 days after explosion). The right panel compares 
the broadband light curves of the model calculated using the Autostructure line data (solid lines) and the Kurucz CD23 linelist (dashed 
lines) . 

(CNdll/CFcIl) 2 ~ 22. 

The variation of the mean opacity with temperature 
(Figure [6]) shows several bumps which reflect changes in 
the ionization state. As the temperature increases, the 
excited levels become more populated, and the number 
of optically thick lines increases. The opacity therefore 
increases with temperature until the gas becomes hot 
enough to ionize. This leads to multiple maxima in the 
mean opacity curve, each of which occur around the tran- 
sition temperatures of the various stages of ionization. At 
sufficiently low temperatures, when the element becomes 
neutral, the opacities cut off sharply, and drop exponen- 
tially with decreasing temperature due to the Boltzmann 
factor in the excited state level populations. 

An important property of the lanthanides is that, rel- 
ative to the iron group, the opacity remains high at rel- 
atively low temperatures. This is because the ioniza- 
tion potentials of the lanthanides are generally ~ 30% 
lower than those of the iron group (see Table 1). For 
neodymium, the mean opacity peaks at T « 5000 K, 
when the ion is mostly singly ionized and cuts offs at 
T < 2500 K when Nd becomes neutral. In comparison, 
the opacity peak for iron occurs at T f=s 7000 K and the 
neutral cutoff is at T < 3500 K. The general persistence 
of the lanthanide opacity to lower temperatures has an 
important impact on the color of the emergent spectra, 
contributing to cooler, redder photospheres. 

Another important feature of the lanthanide opacity is 
the wavelength dependence - while the opacity decreases 
to the red (as there are more lines at bluer wavelengths) , 
the decrease is much slower than that of the iron group 
(Figure [7]). This is due to the much denser energy level 
spacing of the lanthanides, resulting in a much larger 
number of ~ 1 eV optical/infrared transitions. The shal- 
lower opacity profile means that the lanthanides can line 
blanket not only UV wavelengths, but the entire optical 
region of the spectrum. This will influence the color of 
r-process SNe, as photons will eventually be re-emitted 
or fluoresce (through the many lines) to infrared wave- 
lengths where they may escape more easily. 

As seen in Figure [7j the opacity of osmium (Z = 76) is 



10 15 
excitation energy (eV) 

Fig. 5. — Histogram of the number of atomic levels versus 
level energy (bin size = 0.25 eV) in our Autostructure models, 
which illustrates the much greater complexity of the lanthanide 
neodymium (with an open f-shell) as compared to iron (open d- 
shell) and tin (open p-shell). 

very similar to that iron, despite the much higher atomic 
number. This is not surprising, as osmium is a homo- 
logue of iron, with a nearly half open d-shell. Similarly, 
the opacity of the lanthanide cerium (Z = 58) is compa- 
rable to, though slightly less than, that of neodymium. 
This confirms that species with similar complexity mea- 
sures have roughly similar opacities, which we use to de- 
rive approximate opacities for r-process mixtures (j|6]). 

5.1. Uncertainties and Comparison to Existing Data 

Our derived opacities must possess some error, since 
the Autostructure model energies do not exactly match 
the experimental values (Figure[2]) . To estimate how sen- 
sitive the results are to the detailed level energy structure 
and configuration ordering, we examined the Ndll opaci- 
ties derived from the three different optimization schemes 
described in Sj3j The resulting variation provides an es- 
timate of our level of uncertainty. 
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Fig. 6. — Planck mean expansion opacities for three different 
elements, showing the expected dependence on atomic complexity. 
The Nd opacities (blue line, Z = 60, open f-shell) were derived 
from Autostructure models, while the silicon (red line, Z = 14, 
open p-shell) and iron (green line, Z = 26, open d-shell) opacities 
used Kurucz line data. The calculations assume a density p = 
10~ 13 g cm -3 and a time since ejection t e j = 1 days. 
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Fig. 7. — Wavelength dependent line expansion opacities re- 
sulting from Autostructure derived linelists. The opacity of the 
lanthanides (Nd, Ce) is much higher then iron and its d-shell ho- 
mologue, osmium, especially in the infrared. 

Figure [8] shows that the opacities calculated using the 
optl and optS models are quite similar, while the opt2 
model opacities are lower by a factor of ~ 5 at some 
wavelengths. The opt2 model has relatively higher en- 
ergy levels, and hence smaller excited state LTE level 
populations, which is presumably the reason for the lower 
opacities. The optl and opt3 models had similar level en- 
ergies, but the ground state configuration and ordering 
were different. These results suggest that what matters 
most to the opacities is the energy level spacing, and 
not the exact configuration ordering. Given that the 
low lying optS Ndll level energies fairly well reproduce 
experiment, we suspect that further fine tuning of the 
Autostructure model is unlikely to change the result- 
ing opacity by much more than a factor of ~ 2. 

We have also compared our Autostructure opacities 



Fig. 8. — Variations in the wavelength dependent expansion 
opacity for pure neodymium (Z = 60) ejecta obtained using dif- 
ferent Autostructure optimization approaches. These calculations 
adopt a density p = 10 -13 g cm" 3 , temperature T = 4000 K, time 
since ejection t a j = 1 days, and a wavelength binning AA = 0.01A. 



to existing line data from the VALD database, which 



collects ato mic data from a variety of sources (Hcitcr 



et al. 2008). The only high-Z ions with enough lines m 
VALD to derive expansion opacities are Cell and Celll, 
which have wavelengths an d oscillator strengths calcu- 
lated by the Mons group (|Biemont et al 1999 Palmeri 
20001 iQuinet & Biemont 2004 {"The approach of 



et al 2000| |Quinet fc Biemo nt 2004) 
the Mons group to atomic structure is the same as that 
of Kurucz, viz. calculations with Cowan's code utilizing 
extensive experimental energies. In Figure[9j we compare 
the expansion opacities of Ce calculated using the VALD 
linelist and our own Autostructure list. The agreement 
in both the mean and wavelength dependent values is 
good to a factor of ~ 2. Our conclusions about the size 
and wavelength dependence of the lanthanide opacities 
are therefore confirmed when using radiative data from 
independent structure calculations. 

6. OPACITIES OF R-PROCESS MIXTURES 

Although we have only calculated atomic structure 
models for a few ions, the results (Figure [7]) suggest that 
ions of similar complexity have roughly similar opacities. 
This allows us to construct approximate r-process mix- 
tures based on the representative cases. 

In an r-process mixture, the abundance of any indi- 
vidual lanthanide is relatively low (< 1%). Nevertheless, 
these species likely dominate the total opacity. In fact, 
the opacity will depend rather weakly on the exact lan- 
thanide abundance. This is because for the conditions 
found in NSM ejecta, many of the strong lanthanides 
lines are extremely optically thick (r s 1). Such lines 
contribute equally to the expansion opacity regardless 
of the ion's abundance, just as long as that abundance 
remains high enough to keep r s above unity. 

We illustrate this weak dependence on lanthanide 
abundance in Figure [TUJ by computing the opacity of 
a mixture of neodymium and iron. Decreasing the Nd 
mass fraction by a factor of ten (from 100% to 10%) 
only reduces the total opacity of the mixture by ~ 40%. 
Decreasing the Nd mass fraction by two orders of magni- 
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Fig. 9. — Comparison of the expansion opacity for pure cerium (Z = 58) ejecta, computed using the Autostructure line data (red 
lines) and the VALD linelist (green lines). The calculations assume a density p = 10~ 13 g cm -3 , time since ejection t e j = 1 days, and a 
wavelength binning AA = 0.01A. Left: Line expansion opacity versus wavelength for a temperature T = 5000 K. At optical wavelengths, 
the Autostructure results are in reasonably good agreement with the VALD; both are orders of magnitude higher than the opacity of pure 
iron ejecta (brown line). Right: Planck mean opacity as a function of temperature. Only Cell and Celll are included in the calculation. 

this case only Fe and Nd) and the quantity £z specifies 
the total number of elements represented by each. For 
Nd, t;z = 14 to account for all 14 lanthanides, while for 
iron £z = 30 to account for all d-shell elements between 
21 < Z < 80. The quantity pz — Xzp is the density 
of the representative elements, where Xz is the mass 
fraction of each. As an illustrative r-process mixture, we 
assume that the average mass fraction of each lanthanide 
is Xj — 1%, and the average d-shell species fraction is 
Xd = 2%. The remainder of the composition was taken 
to be calcium (s-shell) as a neutral filler. Wc used the 
Nd line data from the optS structure model, and the iron 
line data from the Kurucz CD23 list. 

The dashed line in Figure [10] shows the Planck mean 
opacity of our approximated r-process mixture. Because 
each of the 14 lanthanides is assumed to contribute inde- 
pendently in the sum, the total opacity is essentially 14 
times that of the mixture with only 1% Nd. At certain 
temperatures when the lanthanide opacity dips, the d- 
shell opacity makes a comparable contribution. We note 
that opacity of the mixture can approach the saturation 
level discussed in |2.4[ such that our assumption that the 
strong lines do not overlap can be called into question. 

7. SPECTRA OF NS MERGER EJECTA 

To illustrate the general effect of our r-process opacities 
on the emission from NSM ejecta, we have calculated 
model spectra usin g the Sedona radiation transport code 
( Kasen et al.|2006| . A more comprehensive discussion of 
the light curves and colors of these transients, a nd their 
dependence on the ejecta parameters, is given in [Barnes] 
fe Kasen|p013| . 

As a simple, fiducial ejecta model, we considered a 
spherically symmetric, homologously expanding remnant 
with a broken power-law density profile. The total ejecta 
mass was taken to be M c \ = 0.01 Mm and the ki- 
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Fig. 10. — Dependence of the mean expansion opacity on the 
abundance of lanthanides. The solid lines show the Planck mean 
opacity for various mass fractions of neodymium in a mixture with 
iron. The dashed line shows the opacity of the approximate r- 
process mixture (with all 14 lanthanides) discussed in jj6] 

tudes (from 100% to 1%) reduces the total opacity of the 
mixture by a factor of 5. We find that the Nd opacity 
dominates over that of iron as long as its mass fraction 
is > 10~ 4 . 

The actual r-process ejecta from NSMs will be a het- 
erogeneous mixture of many high Z elements. This mul- 
tiplicity of species should enhance the opacity, as each 
ion contributes a distinct series of lines. To estimate the 
opacity of the mixture we assume the line data of Nd is 
representative of all f-shell species (the lanthanides) and 
that iron is representative of all d-shell elements. We 
ignore the s-shell and p-shell elements since their opac- 
ities will be very low. We then construct the expansion 
opacity of the mixture by generalizing eq. [9] 



pct Ci 



AA, 



1 - exp[-n(pz)}) (15) 



where the sum Z runs over the representative ions (in 



netic energy E = l/2Af e jU^ with a characteristic velocity 
v c = 0.1c. The transport calculations assume the ioniza- 
tion/excitation state is given by LTE, and that the line 
source function is described by the Planck function, i.e., 
the medium is purely absorbing. In reality, the probabil- 
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Fig. 11. — Synthetic spectra (2.5 days after mass ejection) of the 
r-process SN model described in the text, calculated using either 
Kurucz iron group opacities (black line) or our Autostructure de- 
rived r-process opacities (red line). For comparison, we overplot 
blackbody curves of temperature T = 6000 K (black dashed) and 
T = 2500 K (red dashed). The inset shows the corresponding 
bolomctric light curves assuming iron (black) or r-process (red) 
opacities. 

ity of absorption in lines may be small, with fluorescence 
being the more likely result of line interactions. How- 
ever, supernova transport calculations have shown that 
for complex ions, repeated fluorescence among a mul- 
titude of lines has an approximately thermal character 
( [Pinto fc Eastman||2000||Kasen||2006[ ). We use the opac- 
ities ot an r-process mixture derived in £j6l 

The high opacity of r-process material has a signifi- 
cant impact on the predicted radioactive transients from 
NSMs (Figure 111. Compared to previous calculations 
(which assumed iron-like opacities) the predicted bolo- 
metric light curve is of much longer duration, ~ 1 week 
as opposed to ~ 1 day. This is due to the longer effective 
diffusion time through the more opaque ejecta. The peak 
luminosity is also lower, as the radiation suffers greater 
loses due to expansion over this period. The spectrum at 
2.5 days after merger is much redder, with most of the 
flux emitted in the near infrared (~ 1 /mi). Due to the 
extreme line blanketing at bluer wavelengths, the pho- 
tons are eventually redistributed (through lines) to the 
infrared, where the opacities are lower and radiation can 
escape more readily. 

Other than the unusually red color, the r-process spec- 
tra generally resemble those of ordinary SNe, and in 
particular those with high expansion veloc ities (e.g., the 



hyper -energetic Type Ic event, SN 1998bw (Galama et al. 
1998| ). The continuum flux, which is produced by emis- 
sion m the Doppler-broadened forest of lines, resembles a 
blackbody with a few broad (~ 200 A) spectral features. 
It is not easy to associate these features with either ab- 
sorption or emission from a single line; instead they arise 
from blends of many lines. Because our atomic structure 
models do not accurately predict line wavelengths (and 
we only include lines of Nd and Fe) , the location of the 
features in our synthetic spectra are not to be trusted. 
Nevertheless, the general appearance of the model spec- 
tra are likely qualitatively correct. One can anticipate 
where features are most likely to appear by examining the 
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Fig. 12. — Synthetic spectra time series of the r-process SN model 
described in the text. The times since mass ejection are marked 
on the figure. 

energy spacing of the low lying levels of the lanthanides. 

Figure 12 shows the time evolution of the synthetic 
spectra, At the earliest times (< 0.25 days after ejec- 
tion) some flux emerges at optical wavelengths, but this 
phase is short lived. By day 0.5, the optical emission 
has faded, and the spectra evolve relatively slowly there- 
after, with effective blackbody temperatures steady in 
the range T « 2000 - 3000 K. The temporal evolu- 
tion can be understood by considering the mean opac- 
ity curves (e.g., Figure Ejb. At early times, the ejecta is 
relatively hot (> 4000 T\) throughout, and the opacity 
is roughly constant with radius. By day ~ 0.3, however, 
the outermost layers have cooled below < 3000 K, and 
the r-process opacities drop sharply due to lanthanide 
recombination. The ejecta photosphere forms near the 
recombination front (as overlying neutral layers are es- 
sentially transparent) which regulates the effective tem- 
perature to be near the recombination temperature. This 
behavior is similar to the plateau phase of the (hydrogen 
rich) Type IIP SNe, although in this case the opacity 
is due to line blanketing, not electron-scattering. More 
importantly, the temperature at the recombination front 
(Tj ~ 2500 K) is a factor of ~ 2 lower for r-process ejecta, 
as the ionization potentials of the lanthanides (~ 6 eV) 
are lower than that of hydrogen (~ 13.6 eV). It is pos- 
sible that non-thermal ionization from radioactive decay 
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Fig. 13. — Effect of varying the atomic structure model on 
the predicted observables of an r-process SN. The figure shows 
the synthetic spectrum (2.5 days after ejection) calculated using 
Autostructure linelists under different optimization schemes. The 
dashed lines show, for comparison, blackbody curves of tempera- 
ture T = 2500 K. 



products (see { 2.2 ) become important at these phases 
and prevent the outer layers from ever becoming com- 
pletely neutral. 

Our calculated SED's are somewhat sensitive to the 
atomic str ucture model used to generate the line data. 
Figure [13] compares calculations using line data from the 
different Autostructure optimization runs (optl, opt2, 
and optS). The observed differences can be taken as some 
measure of the uncertainty resulting from inaccuracies in 
our atomic structure calculations. Notably, the spectrum 
calculated using the opt2 linelist has significantly higher 
flux in the optical (~ 6000 A) . This is presumably due to 
the lower overall opacity of the opt 2 model (Figure [8]). 
Given the superior match of the optS model to the exper- 
imental level data, we consider the spectral predictions 
using this line data to be the most realistic, however it 
is clear that some significant uncertainties remain. 

Another concern for the spectrum predictions is the 
pote ntial breakdown of the Sobolev approximation (see 
^2. 4). At bluer wavelengths, the mean spacing of strong 



lines can become less than the intrinsic (presumed ther- 
mal) width of the lines, which violates the assumptions 
used to derive an expansion opacity. It is not immedi- 
ately clear how this will impact the results. Two lines 
that overlap exactly will behave like a single line, which 
suggests that an overlapping line should be discounted, 
not double counted. In this case, the opacity should 
saturate at a maximum value K sa t (equation 111. On 



the other hand, in this saturation limit, photons can no 
longer escape lines by redshifting past them, as there are 
no longer any optically thin "windows" between lines. 
The individual line optical depths then become relevant, 
and the effective opacity may be larger than one would 
estimate from the expansion formalism. In practice, the 
impact of line overlap may not be so dramatic - the opac- 
ity generically declines to longer wavelength, and the net 
effect of the transport is to distribute photons to the 
red/infrared where saturation may no longer be an is- 
sue. To fully address the question will require transport 
calculations that dispense with the Sobolev approxima- 



tion and resolve individual line profiles. 

8. SUMMARY AND CONCLUSIONS 

The opacity of r-process ejecta is orders of magnitude 
higher than that of ordinary supernova debris, a fact we 
have demonstrated using new atomic structure calcula- 
tions (Figure [7]), pre-existing line data (Figure [9]), and 
simple physical counting arguments (Figure [TJ. There 
are two physical reasons for the high opacity: 1) Com- 
plexity: The r-process composition includes rare elements 
with complex valence electron structure, in particular the 
lanthanides which have an open valence 4/-shell. Such 
elements have a significantly greater number of levels 
and lines, which results in an overall higher expansion 
opacity; 2) Multiplicity: The r-process produces a het- 
erogenous mixture of many elements, each of which con- 
tributes a distinct series of lines. Since the expansion 
opacity depends on the total number of strong lines (and 
not the strength of any individual line) this diversity of 
the mixture enhances the opacity relative to a more ho- 
mogenous composition. 

While our r-process opacity calculations offer a signif- 
icant improvement over previous estimates (which were 
little more than educated guesses) several uncertainties 
remain. The Autostructure models only approximate 
the level structure of the high-Z elements, and so do 
not correctly predict the wavelength of individual tran- 
sitions. Fortunately, the pseudo-continuum opacities de- 
pend only on the statistical distribution of lines and are 
fairly robust, although our numerical experiments indi- 
cate uncertainties at the factor of ~ 2 level. In the future, 
we can iteratively tune the structure models to better 
reproduce the observed energy levels, although this is a 
time consuming process. Moreover, many of the high-Z 
ions lack good experimental level data. 

Another, more important, uncertainty is that we have 
used the radiative data for one species (Nd, Z — 60) to 
represent all lanthanides. In fact, not all lanthanides are 
equal - the ions whose valence f-shell is nearly open (e.g., 
La, Z = 56) or closed (Yb, Z = 70) are considerable less 
complex, and should have correspondingly lower opaci- 
ties. On the other hand, gadolinium (Gd, Z — 64) has a 
nearly half filled f-shell, and is one of the most complex 
species on the periodic table. Simple counting arguments 
(Figure [I]) suggest that Gd may have an opacity > 10 
times that of Nd. We therefore suspect that our current 
opacities underestimate the true values for an r-process 
mixture where Gd is present at the ~ 1% level. In future 
work, we will quantify the line data for all lanthanides. 
The actinides should also be considered, as they are open 
f-shell as well, although only a few species (e.g., uranium) 
likely have high enough abundance to make a difference. 
In total, the atomic data required for r-process opacities 
is massive, involving numerous structure calculations and 
many billions of lines. 

The high opacity of r-process material has significant 
implications for discovering and interpreting the radioac- 
tively powered transients associated with NSMs. With 
more realistic opacities, the predicted light curves are 
of longer du ration and dimmer at peak (see |Barnes fc 
Kasen||2013[ ). Perhaps even more important to detection 
is that the SED is shifted into the infrared, peaking at 
r~> 1 /im. This is due to the strong line blanketing at 
optical wavelengths, which pushes the photosphere out 
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to cooler layers (T ~ 2500 K) where the lanthanides re- 
combine and become more transparent. Other than the 
very red color, the spectra qualitatively resemble those 
of other high- velocity SNe, with a pseudo blackbody con- 
tinuum and broad (~ 200 A) spectral features. 

The predicted emission at optical wavelengths is some- 
what sensitive to the details of the opacity and its associ- 
ated uncertainties. Spectra calculated using our opt2 ra- 
diative data were fairly bright in the V-band (~ 5000 A) , 
while those calculated using optS data had almost no 
flux at these wavelengths. Because the level structure 
of our optS model agrees better with experiment (and 
given that the overall opacity may be even higher than 
our Nd-based estimates) we consider the latter case to be 
the more likely reality. However, the opacities are still 
not fully converged, and we cannot altogether rule out 
the possibility that r-process SNe may emit some persis- 
tent emission in the optical. 

These results suggest that (to the extent possible) it is 
worthwhile to search for and/or follow up gravitational 
wave sources at red or infrared wavelengths. Optical 
surveys, however, will still be sensitive to a radioactive 
transient if some component of the ejecta is nearly lan- 
thanide free. In fact, it is likely that, in addition to 
the tidal tails, a second component of lighter elements 
(Z < 50) is ejected from a post-merger disk wind. If this 
wind includes radioisotopes with appropriate half-lives 
(e.g., 56 Ni), the light curve may be relativel y bright and 
peak in the optical ( Barnes fc Kasen]|2013[ ). A detailed 
understanding of the wind nucleosynthesis and mixing 
with the tails is important, as our results suggest that 
contamination by lanthanides at just the 10 -3 level may 
significantly raise the opacities. 

Beyond detection, a significant observational challenge 
will be confirmation that a transient is indeed due to a 
NSM. There are likely many classes of stellar explosions 
that produce low mass ejections of radioactive material 



(e.g., Bildsten et al.| 2007 Moriya et al. 2010), and an 



increasing number of fast, faint transients ha ve been ob 
served at optical wavelengths (e.g., [K asliwal et al.||2010l 
Perets et al.|2010||Kasliwal et al.|2012||Eoley et al.|2012n " 
.fortunately, our calculations demonstrate that the opti- 
cal properties of r-process ejecta differ dramatically from 
transients due to lower mass isotopes. A key distinguish- 
ing features is that the r-process SED peaks in the in- 
frared with a nearly constant color temperature regulated 
to the lanthanide recombination temperature, ~ 2500 K. 
The population of brief, infrared variables is mostly un- 
known, but it is possible that the infrared transient sky 
is much cleaner than the optical one. 

By comparing observations of an r-process SN to light 
curve models, one could presumably constrain the mass 



of heavy nuclei ejected in an compact object merger, 
which would go a long way to understanding the un- 
known site(s) of the r-process. One would like to go 
further and spectroscopically study the abundance dis- 
tribution of the outflows. That will be challenging - the 
lines are heavily blended and we do not yet have accurate 
wavelengths for most of them. In the future, though, we 
can refine the line data by tuning structure models to 
match experimental data (where available) and can use 
radiative transport calculations to quantify how global 
abundances variations affect the blended features. While 
measuring a detailed abundance pattern will be hard, 
spectroscopic modeling should permit strong constraints 
on the amount and gross composition of ejecta. 

This assumes that radioactive r-process transients exist 
and that we can find them. Admittedly, we test danger- 
ous waters any time that, lacking observational input, we 
attempt to describe a new astrophysical phenomenon on 
purely theoretical grounds. The situation here is a step 
more treacherous; not only must we rely on simulations 
of a complex macroscopic system, even the microscopic 
structures of our ions are model-based. Obviously, ob- 
servational input is needed to ground the theory. In the 
meantime, the numerical experiment presented in |4.2| 
may offer a bit of comfort. In that example, we calcu- 
lated the light curve of an obliviously crude model of 
a SN la (a uniform 56 Ni blob) using opacities derived 
entirely from ab-initio atomic structure models. Despite 
deliberately ignoring decades of work in the field, our "su- 
pernova from scratch" predictions weren't all that bad, 
and certainly good enough to allow one to search for and 
identify a thermonuclear event. There is then reason to 
hope that the predictions for r-process SNe are (or will 
soon be) good enough for us to know one when we see 
one. 
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